REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


2a.  SECURITY  CLASSIFICATION 


6c.  ADDRESS  (City.  State  and  ZIP  Code) 

College  Park,  Maryland  20904 


3.  DISTRIBUTIQN/AVAILABILITY  OF  REPORT  ,  .  ,  ,  .  „ 

Approved  For  public  release;  distribution 
unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

APOtfftlX*  8  9  -  1  f  f» 


7*.  NAME  OF  MONITORING  ORGANIZATION 

AFOSR,  Dept,  of  Mathematics  and  Information 
Sciences 


7b.  ADORESS  (City,  State  and  ZIP  Code ) 

Bolling  Air  Force  Base 
Washington,  D.C.  20332 


8*.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

AFOSR 


8c.  AOORESS  (City.  State  and  ZIP  Code ) 

Bldg.  410 

Bolling  AFB ,  D.C.  20332 


11.  TITLE  (Include  Security  Classification! 

_ RESEARCH  EQUIPMENT  _ 

1  •'  PERSONAL  AUTHORIS) 

.B.  Kello 


13b.  TIME  COVERED 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable > 

AF0SR-83-0286 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 

element  no. 

61102F 


PROJECT 

NO. 

2304 

TASK 

NO. 

A5 

WORK  UNIT 
NO. 


14.  DATE  OF  REPORT  (Yr,  Mo.,  Da y)  15.  PAGE  COUNT 

TO  1731/55  85/3/25  6 


COSAT  I  COOES 


GROUP  I  SUB.  GR. 


18.  SU8JECT  TERMS  (Continue  on  reverte  if  neceuary  and  identify  by  block  number) 

Numerical  analysis,  finite  elements,  continuation, 
scattering 


19.  ABSTRACT  (Continue  on  reverte  if  neceuary  and  identify  by  block  number ) 

- *  The  purpose  of  the  report  was  to  provide  additional  equipment  for  a  computational 

laboratory  dedicated  to  research  in  numerical  analysis  and  applied  mathematics.  The 
report  summarizes  the  equipment  acquired  and  the  research  being  presently  carried  out  in 
the  laboratory.  The  research  is  on  topics  in  boundary  integral  equations,  adaptive 
finite  elements,  numerical  methods  for  a  scattering  problems,  continuation  methods,  and 
very  accurate  arithmetic  calculations.  /  !  , 

i  -:?•  )  (■ — 


90  01  04  0  46 


31.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


33s.  nami  of  responsible  inoivioual 

Dr.  Robert  N.  Buchal 


33b.  TELEPHONE  NUMBER 
llmeludt  Am  Cod*/ 

(202)767-4939 


33c.  OFFICE  SYMBOL 

NM 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


APOWt  nt*  8  9-1  ?t9 

Final  Report 
Grant  AFOSR  83  0286 
1  September  1983  -  31  January  1985 

Title 

RESEARCH  EQUIPMENT 


Principal  Investigator 


R.  Bruce  Kellogg 

University  of  Maryland 
College  Park,  MD  20742 


o 


rmai  Report  for  Equipment  Grant  AF0SR-83-0286 ,  Sept.  1,1383  -  Jan.  31,  19b5 

The  purpose  of  the  grant  was  to  provide  additional  equipment 
for  a  computational  laboratory  dedicated  to  research  in  numerical 
analysis  and  appMed  mathematics.  The  initial  funding  for  the  laboratory 
was  provided  by  the  National  Science  foundation. 

The  funds  provided  bv  this  grant  were  used  to  acquire  ar,  additional 
computational  node,  ano  additional  memory,  communications  capat i I 1 1 1 es . 
and  graphics  vapat'iites  for  t  n  e  laboratory.  A  detailed  list  <■  t  the 
equipment  aoqmied  by  this  grant  is  as  follows. 
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Following  are  brief  descriptions  of  the  research  being  ■.  aii  ied 
out  in  the  computat  •  one  i  laboratory. 


Project:  Implementation  of  the  spline-trigonometric  method 
Investigator:  Douglas  N.  Arnold  (Department  of  Mathematics) 

A  method  designed  by  the  investigator  for  the  numerical 
solution  of  boundary  integral  equations  is  being  imp  I oriented  and 
tested.  The  method  is  a  Galerkin  method  with  spline  trial 
functions  and  tr  i  gonometr  ic  test  functions.  Among  the  issues 
being  investigated  are  the  proper  numerical  quadrature 
procedures  to  ose  »  th  the  spline-trig  method,  the  use  of  the 
fast  Fourier  transform  m  the  implementation,  the  dependence  of 
the  convergence  on  the  shape  of  the  domain,  and  the  local 
convergence  behaviir.  At  present  a  first  version  of  the  code 
has  been  written. 

This  work  is  partial  I y  supported  by  the  National  Bureau  of 


Standards. 


Project:  Numerical  methods  for  scattering 

Investigators:  R.  B.  Kellogg  UPST) 

A.  K.  Aziz  (IPST) 

Numerical  methods  are  being  investigated  for  the  solution  of 
time  harmonic  scattering  problems  for  moderatelv  high  f requenc 1 es , 

Petweem  the  high  f requeues'  region  where  geometric  opt  in?  techniques 
prevail  and  the  low  frequency  region  where  the  usual  finite  element  and  tin  it 
difference  techniques  are  appropriate.  The  present  stage  of  the 
investigation  involves  the  testing,  on  a  model  problem,  of  some 
modified  finite  element  schemes  that  incorporate  rap'dly  oscillating 
functions  into  the  finite  element  subspace. 

This  result  .s  supported  in  part  bv  the  navy,  through  the  I  Ft  fund  of 


NSWC . 


project:  Path  following  algorithms. 

Investigator:  R.  B.  Kellogg  (IPST) 

Numerical  methods  are  being  developed  and  studied  for  the  solution 
of  n  equations  in  n  +  1  unknowns.  A  path  following  program  has  been 
written,  and  >s  being  tested  on  various  problems  using  the  Apollo 
system.  The  present  stage  of  the  project  involves  using  the  program  to 
follow  curves  generated  by  the  reduced  solution  of  a  rather  complicated 
singular  perturbation  problem  for  a  system  of  differential  equations, 
with  several  boundary  lavers.  The  problem  arises  in  the  modelling  of 
certain  transport  processes  in  physiology. 


Project:  Adaptive  approaches  and  a-posteriori  estimates  for 

solving  differential  equations 

Investigator:  Ivo  Babuska 

The  work  is  partially  supported  by  CNR  Contract  No.  0Q014-85-K-0169 
and  NSF  Grant  DMS-8315316 

1)  The  p  and  h-p  versions  of  the  finite  element  method. 

The  theory  of  the  p  and  h-p  versions  of  the  finite  element  method 
in  one  and  two  dimensions  was  developed.  The  practical  applicability  of  the 
theoretical  conclusions  was  tested  by  the  numerical  treatment  of  relatively 
large  scale  problems  of  structural  mechanics.  The  results  have  been  already 
used  in  the  commerical  p-version  finite  element  program  PROBE  (Noetic 
Technologies,  Inc.).  The  numerical  results  are  part  of  two  Ph.D.  theses 
awarded  in  1985. 

2)  Numerical  solution  of  linear  and  nonlinear  boundary  value  problems 

for  ordinary  differential  equations. 

The  theory  of  the  factorization  method  was  developed  and  a  general 
program  was  written.  The  program  was  tested  on  various  problems  of  an 
engineering  character  as,  for  example,  the  problem  of  the  stability  of  a 
towed  underwater  cable.  The  results  are  related  to  the  Ph.D.  theses  awarded 
in  1984. 

3)  Adaptive  finite  element  solution  for  nonlinear  problems. 

The  previous  work,  which  resulted  in  the  design  and  implementation 
of  the  adaptive  solver  FEARS  for  linear  partial  differential  equations,  is 
being  extended  to  the  design  of  the  nonlinear  solver  NFEARS. 


Project:  Level-index  arithmetic 

Investigators:  F.  W.  J.  Olver  (Research  Professor,  IPST)  and 
D.  W.  Lozier  (Adjunct  Professor,  IPST). 

The  level-index  (Zl)  number  system  is  being  investigated  as  an  alternative 

to  the  floating-point  system  for  internal  use  on  computers.  In  this  system  a 

nonnegative  real  number  x  is  represented  as  its  generalized  logarithm 

Z  +  l n'  'x  ,  where  Zn(  'x  is  the  Z  th  iterated  logarithm  of  x  ,  the  integer 

(Z) 

t  being  determined  by  the  condition  0  £  Zny  'x  <  1  .  The  main  advantage 
of  the  Zi  system  would  be  to  eliminate  overflow  and  underflow  problems.  The 
Apollo  will  be  used  to  simulate  the  Zl  system  in  software  in  order  to  test  the 
performance  of  the  system  on  various  standard  algorithms  in  numerical  analysis, 
particularly  ones  that  are  known  to  be  prone  to  overflow  or  underflow.  A 
compiler  is  nearing  completion. 

This  work  is  supported  in  part  by  the  V.  S.  Army  Research  Office,  Durham 
under  Contract  DAAG  29-84-K-0022. 

Reference:  C.  W.  Clenshaw  and  F.  W.  J.  Olver  "Beyond  floating  point",  J.  A.C.M. 
31,  319-328  (1984). 


